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Abstract 



This paper considers the problem of spectrum sensing in cognitive radio networks when the primary 

O . user employs Orthogonal Frequency Division Multiplexing (OFDM). We develop cooperative sequential 

detection algorithms based on energy detectors and the autocorrelation property of cyclic prefix (CP) 

K^ I used in OFDM systems and compare their performances. We show that sequential detection provides 

^O ■ much better performance than the traditional fixed sample size (snapshot) based detectors. We also study 

m 

the effect of model uncertainties such as timing and frequency offset, IQ-imbalance and uncertainty in 
noise and transmit power on the performance of the detectors. We modify the detectors to mitigate the 
effects of these impairments. The performance of the proposed algorithms are studied via simulations. 
It is shown that energy detector performs significantly better than the CP-based detector, except in case 



^^ i of a snapshot detector with noise power uncertainty. Also, unlike for the CP-based detector, most of the 

C^ ' above mentioned impairments have no effect on the energy detector. || 
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I. Introduction 

Cognitive Radios, also called secondary users, use the radio spectrum licensed to other (primary) users. 
They perform radio environment analysis, identify the spectral holes and then operate in those holes [1]. 
Spectrum Sensing (SS), an essential first step in enabling Cognitive Radio (CR) technology, involves 

(1) identifying spectrum opportunities by detecting holes (white space) when they become available and 

(2) detecting when the primary reclaims an identified spectral hole. This needs to be done such that the 
guaranteed maximum interference levels to the primary are not exceeded and there is efficient use of 
spectrum by the secondary [2]. Thus, the cognitive radio needs to detect reliably, quickly and robustly, 
possibly weak primary user signals. For example, the IEEE 802.22 standard [1] requires a sensitivity 
of -116dBm, while keeping the probability of miss detection under 0.1, using a sensing duration < 2 
seconds. 

Orthogonal Frequency Division Multiplexing (OFDM) is used in 802.1 la/g wireless LAN's (WLAN), 
Wireless MAN's (IEEE 802. 16 WiMAX,) 3GPP Long Term Evolution (LTE), etc. Because of its widespread 
acceptance and deployment, it is likely that a primary user would be using OFDM, thus making the 
problem of detecting OFDM signals especially relevant for cognitive radio. Most of the OFDM systems 
also employ a cyclic prefix (CP), which implies that the autocorrelation is non-zero at delays of the 
useful symbol length — a property that can be exploited for spectrum sensing [3]. Literature on spectrum 
sensing is vast, despite it being a relatively recent topic of research. We now briefly summarize the 
relevant recent work. 

A. Literature Survey 

For spectrum sensing, primarily three signal processing techniques [4] are proposed in literature: 
matched filter [5], energy detection [5] and cyclo-stationary feature detection [6]. Matched filtering is 
optimal but requires detailed knowledge of the primary signal. When no such knowledge is available, 
an energy detector is optimal [5]. Hence, most of the literature is based on energy detection. However, 
unlike for the matched-filter and the cyclo-stationary detectors, it suffers from the so-called SNR wall 
problem in the presence of transmit power or receiver noise power uncertainties ([5], [7]). 

Another important problem encountered in cognitive radios is the hidden node problem caused due to 
shadowing or time-varying multipath fading. To alleviate this problem, cooperative sequential spectrum 
sensing algorithms are suggested. Cooperative spectrum sensing where the decisions of different secon- 



daries are fused to obtain the final decision has been studied in [4], [8], [9] and [10]. Sequential detection 
techniques have been used in [11], [12], [13] and [14]. 

Spectrum sensing in an OFDM environment has been studied in [3], [15] and for Orthogonal Frequency 
Division Multiple Access (OFDM A) systems in [16], [17]. In [15] CP correlation based snapshot and 
sequential detectors are studied. In [16] the effect of time and frequency offset on CP based snapshot 
detector is studied. In [17], the joint sensing and channel scheduling problem is addressed for OFDM- 
based CR systems. 

B. Our Contribution 

In this paper, we provide spectrum sensing algorithms (for detecting spectral holes in time) when the 
primary is using OFDM. As in recent work on spectrum sensing in OFDM, we exploit the autocorrelation 
property in our spectrum sensing algorithms. We compare this with energy detector based algorithms. 
Furthermore, we study how some of the common impairments in a secondary receiver [18] like frequency 
offset, timing offset, IQ-imbalance [19] and uncertainty in the noise power and unknown channel gains 
affect the statistics, and thus the performance of these detectors. We propose techniques to modify the 
detectors to mitigate some of the losses. For simplicity, these techniques are first presented in the context 
of fixed sample-size detectors, but are then extended to the sequential detection setup as well. Here, 
we are primarily interested in the sequential detection algorithms which are more efficient than fixed 
sample size (snapshot) detectors, commonly used in the literature. In particular, we use DualCUSUM, 
a distributed, cooperative sequential-detection algorithm developed in [20]. It was shown in [14] that 
DualCUSUM outperforms several other existing spectrum sensing algorithms. However in [14] it was 
not studied in an OFDM context and only the energy detector was studied. 

We show that, unlike for the CP-based detector, the energy detector is inherently robust to timing 
offset, frequency offset and IQ imbalance. However, it is sensitive to noise power uncertainties, while 
the CP based detector is not. In both the cases, we modify the detector and improve the performance in 
the presence of impairments. Our study shows that energy detector substantially outperforms CP-based 
detectors in the setup considered here (except in case of snapshot detector with noise power uncertainties). 

This paper is organized as follows. Section II describes the OFDM model. We also present our 
cooperative sequential detection based setup. In Section III we study the effect of the impairments on the 
snapshot CP-detector and present possible techniques to overcome the same. In Section IV we consider 



the energy detector in a snapshot setup and study the effect of the impairments. In Section V we extend 
these techniques to the sequential change detection algorithms. Section VI concludes the paper. 

II. Model 

We consider a CR network which is sensing a primary using an OFDM system. The OFDM system of 
the primary (Figure 1) consists of L^ narrowband signals Dq, Di, ..., -Dl^-i carried by the subcarriers. An 
OFDM symbol is obtained by passing the L^ signals through an inverse fast Fourier transform (IFFT). 
In addition, a cyclic prefix of length Lc is also appended to make the total OFDM symbol duration 
Ls = Ld + Lc- The inter-carrier spacing is denoted by A/. 

Do Di DLd-1 



CP 



IFFT 



Fig. 1. OFDM symbol construction. 



Define, d{k) = En=o^ Dnc' ^^-i . The baseband OFDM signal at time k is, 

I d{k + Ld-L^), k = 0,l,-,Lc-l, 

S{k) = <f (1) 

[ d{k-Lc), k = Lc,Lc + l,...,Ls -1. 

By the Central Limit Theorem (CLT), ^(A;) is approximately Gaussian, since it is a linear combination 
of Ld random signals. Also, as E[Dn] = 0, E[S{k)] = 0. We will take Ld = 64, Lc = 16 and 
A/ = ^MHz in our simulations, although our analysis and algorithms are general. These parameters 
are assumed to be known to the secondary systems. 

In ([T]), we observe that the symbols S{0), ..., S{Lc — 1) are repeated as S{Ls — Lc), ■■■, S{Ls — 1). 
Thus, if we correlate this sequence with a shift of Ld, we will get a good correlation in case the primary 
is transmitting; otherwise not. CP based detectors exploit this property in detecting the primary signal 
[3]. 

We consider a cognitive radio system with L secondary users that sense a channel via CP-detectors. 
Later on, we will also consider energy detectors. The observations made on the channel by these 
secondaries are processed and sent to a fusion center, which makes a decision on whether the channel 



is free or not. Then, that decision is sent to all the secondary users for possible use of the channel. 
The secondary system has to detect when the primary starts transmission (OFF— s-ON) and when it stops 
(ON^-OFF). In the following, we explain our setup for OFF— J-ON, but our algorithms work for ON^-OFF 
also. 

Let the primary start transmission at a random time T. Then at time k the signal received by the l^^ 
secondary is, 

Xik.l) = \ ""•■"■ ^=1.2.....T-1, ^^^ 

\ h{l)S{k) + N{k,l), k = T,T + l,.... 

where h{l) is the channel gain of the Z*^ user, and N{k, I) is observation noise at the l*^ user. We 
assume that the fading is frequency flat and remains constant during the interval of observation (say, 
approximately, for a duration of ON/OFF period). Slow fading scenarios with primary staying ON and 
OFF for a few seconds will approximately satisfy this. This assumption is commonly made in the literature 
[9], [12]. (We will see below that most of our algorithms can be used for the frequency selective fading as 
well). We also assume that {d{k), k > 1} and {N{k), A; > 1} are independent and identically distributed 
(i.i.d.) sequences independent of each other and T. Thus, pre-change X{k,l) are i.i.d. with distribution 
Nc{0,a'^ i), where Nc denotes circularly symmetric complex-normal distribution and a^ ; the noise power 
at node / ( N denotes normal distribution). The post-change distribution of X{k, I) is Nc{0, c^ ; + o's i) 
where a^^ is the received power of primary node at node I. The effect of different channel gains is 
absorbed into a^^. 

The aim is to detect the change (at random time T) at the fusion center as soon as possible at a time 
T (> T) (i.e., to minimize E[{t — T)+], where (x)^ = max(0,x)) using the messages transmitted from 
the L secondaries, with an upper bound on probability of false alarm, PpA — P{t < T) < a. For 
this, each of the L nodes uses its observation X{k,l) to generate a signal Y{k,l) and transmits to the 
fusion center. The data received at the fusion center is corrupted by the additive white Gaussian noise 
(AWGN) at the receiver. The fusion center uses the observations Y{k, 1), ... , Y{k, L) to decide between 
the two hypotheses Hq (primary not transmitting) and Hi. If Hq is chosen, the secondaries continue to 
use the channel in slot k and the spectrum sensing session continues (they may use part of each slot 
for sensing and the rest for transmission). If Hi is detected, the secondaries typically switch over to 
an alternate channel. Our algorithms do not change (although the parameters can) if we interchange the 
role of Hq and Hi. To transmit Y{k, 1), ...,Y{k,L) from the L secondaries to their fusion node, they 
need a Multiple Access Control (MAC) protocol. Time Division Multiple Access (TDMA) is the most 



commonly used protocol. 

We have developed (see [11], [20]) a robust cooperative algorithm for spectrum sensing in this setup. 
In this paper, we study this algorithm in the OFDM setup. We also modify it to take care of the various 
impairments commonly encountered in OFDM. 

The rest of the paper is organized as follows. First, in Section III, we study the CP-detector under 
the classical snapshot setup in the presence of different impairments. In Section IV we study the energy 
detector under different impairments. In Section V we adopt the different estimation schemes used in 
Section III and Section IV, to the cooperative sequential detection setup of DualCUSUM to detect the 
presence of the primary signal. Further, we compare the performance of the energy detector with the 
CP-detector. 

III. Cyclic Prefix Based Detector 

In this section, we explain the CP correlation based snapshot detector in the context of a single 
secondary node (thus the subscript I will be omitted in the notation) and present how we mitigate the 
effects of different impairments and uncertainties. Given a number of observations X{1), ... ,X{MLs) 
from M slots of OFDM symbols, we want to detect if Hq or Hi is true. We use the Neymon-Pearson 
(NP) method for detection. We compute the autocorrelation at lag Ld, 

M-l L, 

^=J^J2I1 ^(-^'^^ + ^)^* (-^'^^ + i + Ld) (3) 

" j=0 i=l 

where X* is the complex conjugate of X. The above detector assumes perfect OFDM symbol level 
synchronization and thus correlates only the exact set of samples which would be repeated in the CP 
under Hi. Using the CLT, it can be shown [3] that R ~ A'^c(0, o-g) under Hq and R -^ Nc{(t1, af) under 
Hi, where ag = -^ and af = -^^^^x^- ^^ 1°^ '^^^ (i-^- ^s < ^^)' ^o ^ ^i- We work under this 
assumption, as in CR our main concern is signal detection at low SNR. 

Since the post-change mean is real, under low SNR conditions, detection is based on the real part 

4 

Rr = Real(i?) as Rr ~ N{0,a^) under Hq and Rr ~ iV(o-^,cr^) under Hi, where a^ ^ f^ . The 
detection rule is of the form R^ > X for declaration of Hi. 

In case of frequency selective fading the distribution of R under Hq remains as above. Under Hi, i? is 
still approximately Gaussian with mean and variance now depending on the fading parameters. Knowing 
the fading parameters, one can obtain the detection rule via NP lemma. 

Next, we discuss different impairments and possible techniques to mitigate their effects. 



A. Timing Offset 

Timing offset occurs because the cognitive receiver may not know wliere the OFDM symbol boundary 
starts in the received set of samples. Thus, it may not know the exact set of samples to correlate in ^. 
If it correlates at an incorrect position, E[R] « under Hi. One possible way to take care of this is to 
correlate for the duration of the entire OFDM symbol: 



— Y^X{i)X*{z + L,)\. (4) 



Now a^ sa — - — under either hypothesis. The post change mean under Hi is fia"^, where fi = —. 
Because of this, one can expect the performance to degrade. Using M = 100 and SNR = —WdB 
{a'^ = 20, a'^ = 2), we simulated this setup to show the effects of timing offset. We use these parameters 
throughout this section. The unknown timing offset is chosen as 30 samples. For different values of the 
probability of false alarm pfa (detecting Hi while Hq is true), the detection probability pd (detecting 
Hi while Hi is true) is shown in Table H To regain some of the lost performance, instead of correlating 
over the entire set of samples, we can estimate the unknown timing offset 6* by a maximum likelihood 
estimator (MLE) [18] as: 

9ml = arg max {Re{R{9)) - ujP{9)}, (5) 

6»e{l,2,...,L,-l} 

where 

M-l L, 

R{9) =^Y1 ^(-^'^^ + ' + ^)^*(J^s + i + Ld + 9), 

j=0 i=l 



M-l L 



^(^) = \Y1 flil^i^Ls + i + 9)\^ + \XiJL, + i + La + 9)\^),and oj 



a 



2 



Under low SNR, uj is small. Also, for a large number of OFDM symbols, P{9) Ki MLc{cr'^ + a^^) 
under Hi and P{9) ^ MLcCr^ under Hq. Thus, ujP{9) does not affect the max operation in ([5]), and we 
use the simplified estimator 

9ml = arg max \Real(R(9))} . (6) 

6»g{l,2,...,L,-l} 

This estimator has the advantage of not requiring knowledge of a'^ or o"^. Then we use the decision 
statistic 

/ Af-l L, 

Rr = Real [ml^YI ^(^'^^ + ' + ^ml)X* {jLs + i + L^ + 9 ml) \ (7) 

\ i=o i=i 



Pfa 


No Impairments 

m 


Correlating over 
entire symbol ^ 


Timing Offset 
estimate (|7]i 


0.05 


0.9999 


0.7921 


0.9975 


0.025 


0.9996 


0.7010 


0.9944 


0.01 


0.9988 


0.5770 


0.9880 



TABLE I 

SNAPSHOT CP DETECTOR: EFFECT OF TIMING OFFSET IN pd.THE UNKNOWN OFFSET IS SET AS 30. 

instead of ([3]). Under Hq and Hi, now Rr is no longer normally distributed, but the Gaussian distribution 
still provides a good fit: the empirical distribution of Rr under either hypothesis and the normal tit is 
shown in Fig|2l Thus, we use this approximation for designing the detection threshold and performance 
analysis. However, as the variances under Hq and Hi are different, the optimal likelihood ratio is not a 
linear function of Rr and involves knowledge of a^ at the CR, which is not desirable. Thus, we propose 
to continue to use a test of the form R^ > A which is sub-optimal in this case, and could be viewed as 
a non-parametric test. The performance comparison is shown in Table Jl It can be seen that we recover 
most of the performance lost due to timing offset. 




Fig. 2. Empirical CDFs of Rr under Ho (signal absent) and Hi (signal present) with timing offset. Note that the Gaussian 
approximation provides a good fit to the empirical CDF. 

B. Frequency Offset 

Let us now consider the scenario when only a frequency offset is present (i.e., the timing offset is 
assumed to be known). Let the frequency offset (between the cognitive receiver oscillator and the primary 
transmitter oscillator) be denoted by (/), normaUzed with respect to the carrier spacing A/ . The received 
signal can be written as X{k) = S{k)e^'^^^^ +N{k). Under iJi, R ~ Nc{a'^e~^'^'^'l',aj). If the receiver is 
not aware of the frequency offset, the post change Rr ~ Nc{agCos{2TT(j)) , a^), degrading the performance 



(see Table HIl for cj) = 0.1). To mitigate this effect, we estimate the frequency offset (p via an MLE 0ml- 
The log likelihood ratio can be shown to be proportional to 

2(jf (iir- cos(27r(/)) + Ri sin(27r0))/cj? (8) 

where Rr and Ri are the real and imaginary parts of R, respectively. It can be shown that (/}ml = 
—ZR/2tt, and we use this estimate in the NP test. Thus, the optimal test becomes \R\'^ > A'. Under 
Hq, |i?p has an exponential distribution, and under Hi, it has a non-central Chi-square distribution. The 
performance is shown in Table |lll Note that once again, most of the performance loss is recovered. 
When both timing and frequency offset are present, one can estimate these as 

9ml = argmax|i2(6')|, 4>ml = ---/^R{0ml)- (9) 

9 zvr 

We will use these estimates when we consider all impairments together. 

C. IQ-Imbalance 

IQ-imbalance occurs due to non-ideal front end components in the receiver [19] resulting in the 
amplitude and phase imbalance in the inphase (I) and quadrature (Q) components of the signal. In 
the presence of IQ-imbalance the actual received signal is written as 

X{k) = aY{k) + ^Y*{k) (10) 

where and ^(A;) = S{k) + N{k),a = cos(A(/)) + jesin(A(/));/3 = ecos(A</)) — jsin(A</)) and e and A(/) 
are the amplitude and phase imbalance parameters respectively. It can be shown that in the presence of 
IQ-imbalance, 



A^(0, (t|q), underHo, 



Rr^ { " ' '^" " (11) 

^jn), under Hi 

where ajg ^ o-;j^((l + e^)2 + 4|Cip)/2MLc under low SNR conditions and Ci = a/3*. The performance 
of the detector is shown in Table HIl for A0 = 10°; e = 0.2. We see that the performance of the detector 
degrades slightly even when knowledge of imbalance parameters are assumed but not compensated for. 
However, we can improve performance by compensating for the imbalance. We use the algorithm in [19] 
to compensate for IQ-Imbalance before starting the CP-detector. The imbalance parameters are estimated 
and corrected for as follows. Let 

^2AE!!fllM^ and i^^^. 
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Pfa 


Frequency Offset 
without compensation 


Frequency Offset 
with compensation ^ 


IQ-Imbalance, 
No compensation dill) 


IQ-Imbalance 


with compensation (|12|13|l 


0.05 


0.9965 


0.9989 


0.9991 


0.9999 


0.025 


0.9913 


0.9975 


0.9977 


0.9996 


0.01 


0.9794 


0.9939 


0.9937 


0.9988 



TABLE II 

Snapshot CP detector: pd under Frequency Offset and IQ-Imbalance. The normalized frequency offset 

WAS SET TO 0. 1 AND THE IQ IMBALANCE PARAMETERS A0 AND e ARE SET TO 10° AND 0.2 RESPECTIVELY. 



where Xr and Xi are the real and imaginary parts of X respectively. Then, one can correct the amplitude 
imbalance by 



Zr{k) 



Xr{k) 



Zi{k) 



Xi{k) 



1 + e' "'""' 1-e 

Assuming the phase imbalance € [— 7r/4, 7r/4], it is estimated and corrected as, 



(12) 



-.ML 



Zr"'Xr(i)X,ii) 



A^ 



sm-\2S) 



(13) 



Zf''^{X?ii)+Xf{^)y ' 2 

Then, instead of using the observations X{k), we use X'{k) with real and imaginary components 

X^{k) cos(A0) sin(A0) Zr{k) 

X[{k) sin(A0) cos(A(^) Zi{k) 

for the CP detector. The performance of the detector with this estimator is shown in Table Ull We see 
almost no performance loss. 

From these results, we see that the performance loss due to the IQ imbalance could be ignored. 
However, we have found that it does cause non-negligible degradation when there are other impairments 
mentioned above. Then the improvement resulting from the compensation procedure described by (fT2l ). 
(IT3l) can be more significant. 



D. Noise/Transmit Power Uncertainty 

In a cognitive radio setting, the receiver noise power o"^ and the received signal power a1, may often 
not be precisely known to the CR [5]. We now address the detection problem under these uncertainties. 
Since the variance of Rj. is dependent on the noise power, the detection threshold cannot be set without 
its knowledge at the CR receiver. Thus, the noise power is estimated as 



•.ML 






var{X) 



T^tJrX{^)X^^ 



ML, 



(14) 
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Pfa 


Noise power 
estimation (Il4ll 


All impairments 
(^ with noise power estimation) 


All impairments 
with all compensation 


0.05 


0.9999 


0.5122 


0.9712 


0.02 


0.9995 


0.3908 


0.9556 


0.01 


0.9976 


0.2614 


0.9300 



TABLE III 

SNAPSHOT CP DETECTOR: pd UNDER NOISE UNCERTAINTY AND ALL IMPAIRMENTS 

and this is used to set the threshold A to achieve desired pja- This causes a minor performance loss 
if this estimate is obtained when Hi is true, since, then var{X) ^ (a^ + a'^)'^ /MLg. However at low 
SNR's this causes small estimation error. This can be verified from Table |llll Also, we have been using 
tests of the form Rr > X or \Rr\ > X (partly motivated by the constraints of the present section), and 
the statistics of Rr do not depend upon a"^ under Hq. Thus, knowledge of receive signal power is not 
necessary to set the threshold A to achieve the desired pja- 



E. All Impairments 

In this section, we simulate the performance of the fixed sample size CP-detector when all impairments 
are present. First, the detector estimates and compensates for IQ imbalance using (fT2l) and (fT3l ). Then, 
the variance of received signal is estimated to set the threshold. Next, the optimal timing and frequency 
offsets are estimated using Q and the test is of the form \R\ > A. The performance is shown in Table 
Unl under the impairments and data statistic given in this section. We see that the estimation schemes 
recover most of the losses. 

For reference, we also compare with the detector in (|4]) in the presence of IQ-imbalance, frequency 
offset and noise uncertainty. Noise uncertainty for this detector is taken care of as in Section III.D (i.e., 
estimating the noise variance to adjust the threshold) as this is necessary to set the threshold. It takes 
care of timing offset by correlating over the entire OFDM symbol duration, but the detector is unaware 
of frequency offset and IQ imbalance. Thus, even with partially compensating for the impairments, the 
performance can be very poor. However, from the last column in Table |llll we see that using the methods 
presented here, most of the losses can be recovered. Motivated by these estimation schemes, we mitigate 
the effects of these impairments in the sequential detection algorithm, DualCUSUM, in Section V for 
the CP-detector. 
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IV. Energy Detector 

In this section, we study the performance of the energy detector under a snapshot setup for a single 
secondary node (as in Section EI). We study the effect of different impairments and explore possible 
techniques to mitigate the same. We compute the energy 



ML, 



^ = T7=^>:i^»i'- (15) 



1=1 



Using the CLT, it can be shown that 



y ~ < 



(16) 

N{K + <^'),^^l^ + ^m), under i?i. 



The additional term in variance of V under Hi arises due to the presence of the cyclic prefix. But at 
low SNR assumptions, it is easy to see that V is approximately distributed as r^ N{a1 + o"^, "jy^^™ ) 
under Hi. We work under this assumption. All the likelihood ratio tests are of the form y > A, as the 
likelihood ratio test will involve the knowledge of primary signal power al- 

For the frequency selective case, as in CP detector, V will again be approximately Gaussian with the 
mean and variance under Hi different from the frequency flat case. 

A. Timing Ojfset 

The effect of timing offset in the context of energy detection is that in a set of MLs samples, we do 
not know exactly how many samples would belong to the cyclic prefix portion of the OFDM symbol. 
This in turn implies that we would not know exactly how many of the terms in the expression for V 
given by (fT6l ) would be correlated. For example, for a timing offset of 6* S {L^, ■■■,Ls — Lc — 1}, 

i.e., the second term in variance could be different from that given by under Hi could be different from 
that given by dlTl) . depending upon the value of 6. But under low SNR conditions the effect of this is 
negligible, and thus timing offset does not affect the performance of the energy detector. The results are 
shown in Table |lVl The parameters are M = 40 , SNR = -lOdB (cr^ = 20, a^ = 2) and the unknown 
timing offset was chosen as 30. The number of OFDM symbols used in this section is different from 
that in Section III. This is because, with M = 40 OFDM symbols, the energy detector provides a much 
superior performance compared to the CP detector, under no noise uncertainty. 
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Pfa 


No Impairments 
113J 


Timing offset 


Frequency Offset 


0.05 


0.9999 


0.9999 


0.9999 


0.025 


0.9996 


0.9997 


0.9996 


0.01 


0.9988 


0.9989 


0.9988 



TABLE IV 

Snapshot Energy Detector: pd under timing oefset=30 and normalized frequency offset=0. 1 

B. Frequency Offset 

As the effect of frequency offset is a rotation of X{k) and since the distribution of X{k) is rotationally 
invariant, the statistics of V is not affected by the frequency offset. Thus, the performance of the Energy 
detector is not affected by frequency offset. The assumption here is that the loss in signal energy due to 
the implicit band pass filtering prior to energy detection is negligible. The results are shown in Table JV] 
for a frequency offset of = 0.1 (normalized by the inter carrier spacing A/). 



C. IQ-Imbalance 

In the presence of IQ-imbalance, the statistics of energy detector are: 



V 



(18) 



iV^(l + a(a^ + a,^),i^-^^^lAii-£L-Si^J^J, under //i. 

The performance of IQ imbalance under no compensation and with the compensation scheme of Section 
III.C is shown in Table JV] 



D. Noise/Transmit Power uncertainty 

It is well known that under presence of noise uncertainty, the energy detector has a SNR wall and 
the performance suffers. This is illustrated in this subsection. Let cr^ G [^,(j^(5], where 5 denotes the 
uncertainty level and a^ denotes the nominal noise power used in other sections. The performance of the 
energy detector when a"^ = 10 and 5 = 1.08 (corresponding to 0.33 dB) is shown in Table 5. The energy 
detector sets the threshold for ct^(5 and thus the probability of detection significantly degrades. Also, as 
mentioned in Section III, since the tests are of the form V > \, knowledge of a1 is not necessary to 
meet a desired pja- 
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Pfa 


IQ Imbalance 
(No Compensation) 


IQ Imbalance 
(Compensation) 


Noise Uncertainty 


Timing, Frequency 

and IQ Imbalance 

with compensation for IQ 


All Impairments 

with compensation 

forlQ 


0.05 


0.9992 


0.9996 


0.2785 


0.9995 


0.2563 


0.025 


0.9981 


0.9989 


0.1814 


0.9985 


0.1675 


0.01 


0.9941 


0.9968 


0.1007 


0.9965 


0.0921 



TABLE V 
Snapshot Energy Detector: pd under IQ-Imbalance, noise uncertainty and all impairments for the 

PARAMETERS OF SEC. IV 



Pfa 


CP-no Impairments 
(No Compensation) 


CP-timing, freq offsets 

and IQ Imbalance with 

compensation for all of Sec. Ill 


All Impairments 

including noise uncertainty 

with compensation 

for all of Sec. Ill 


0.05 


0.9606 


0.7097 


0.5701 


0.025 


0.9606 


0.6173 


0.4680 


0.01 


0.8728 


0.4905 


0.3334 



TABLE VI 

Snapshot CP detector: pa under the parameters of Sec. IV 



E. All Impairments 

In this section, we simulate the performance when all impairments excluding noise uncertainty are 
present (IQ-imbalance is compensated) and then later include the effect of noise uncertainty. Also we 
have simulated the performance of the CP-detector for the parameters of this section (i.e. M = 40) in 
Table |Vll Comparing Table Ivl and Table IVll we can see under all impairments excluding noise uncertainty, 
the energy detector has a better performance than the CP detector compare column (iv) of Table |V] and 
column (i) of Table I VII When noise uncertainty is present, the performance of the energy detector degrades 
significantly compared to cyclic prefix detector and thus, in a snapshot setup, the CP-detector is more 
robust to these impairments (last columns of Tables |V] and IVll ) than the energy detector. 
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V. Cooperative Sequential Sensing of OFDM 

The advantages of spectrum sensing by cooperative means, i.e., using multiple nodes to sense the 
spectrum, are well known [4], [10]. Furthermore, sequential detection is also known to perform better 
than snapshot detection. In this section, we apply cooperative sequential detection algorithms developed 
in [11], [14], [20] for sensing of the OFDM signal in the setup of Section II. Interested readers are 
referred to [11], [14], [20] for a more detailed introduction to sequential detection and its advantages. 

We compare the performance of cooperative algorithms with different levels of impairments. Dual- 
CUSUM uses the well known CUSUM algorithm [21] at the cognitive receivers as well as at the fusion 
node for detection of change (ON -^ OFF and OFF -^ ON of the primary). CUSUM is known to be 
optimal in different scenarios and uses the log likelihood ratio. Consequently, DualCUSUM has also been 
shown to perform very well ([14], [20]). In the following, we use DualCUSUM in our present scenario 
and treat both energy detector and cyclic -prefix based detector simultaneously. We use the estimation 
schemes (wherever applicable) discussed in Section III and Section IV (suitably modified), overcoming 
the effects of different impairments. 

In Tables IVIII and IVIIII we provide the performance of DualCUSUM and its variants. The parameters 
used for simulations are described in Section V.E which also compares the algorithms in different 
scenarios. 

A. Dual CUSUM with No Impairments 

This is the ideal scenario where there are none of the impairments mentioned in Section III. For the 
cyclic prefix detector, correlation is done only over the length of samples corresponding to the cyclic 
prefix. Since all the parameters, including noise variance and received primary power are known, one 
can apply the DualCUSUM [20] as explained briefly below. 

1) Each node I computes the log likelihood ratio (LLR) ^j ^ of Rr{j,l) in each slot j{> 1) of Lg 
samples as 

Rrij, I) = Real i -i ^ X{{j - 1)L, + i, l)X*{{j - 1)L, + L^ + i, /) i , (19) 



L,{2aliRr{j,l)-ali) 
2ai 



(j,i = '' : ; — (20) 



and computes the cumulative summation (CUSUM) 

Wj^i = {Wj.i^i + ^j,i)+, Wo,i = 0. (21) 
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2) If the CUSUM crosses a threshold 7, it transmits a message Yj^i = bl^^/..^^^ to the fusion node 
(i.e., it sends a 1 with ampUtude b). 

3) The fusion center receives Yj in slot j where 

Yj=Y.^U + Zr (22) 

I 

and {Zj} is i.i.d. receiver noise with distribution N{0,af,j). 

4) The fusion node also runs CUSUM based on its input Yj by using the log likelihood ratio rjj as 
follows: 

Ffc = (Ffc_i+7?,)+, Fo = 0, rj, = '^^'^i ~ ^^^^ , (23) 

where / is a design parameter. 

5) Fusion node finally declares change at time r if F^ exceeds a threshold j3, i.e., 

r = inf{A: : F^ > f3}. (24) 

The parameters 7, /3, 6, I affect the performance of the algorithm and the techniques developed in [20] 
can be used to optimize performance. One computes EDD = E[{t — T)+] subject to the probability of 
false alarm PpA <a = P[t <T]. 

For the energy detector, the algorithm is the same as the above with minor modifications. The energy 
is computed as 

'''■'''^ = - ML, <«' 

and ^j^i is the LLR computed with pre and post change distributions being N{a'^, a^/MLg) and Nla"^ + 
al, {al + a'if/MLs) respectively, 

^^'^ 2 °^ V {<y1 + ^^)^ ) </ML, {al + alY/ML, ' ^ ^^ 

For frequency selective fading, V{j, I) in (l25l ) will not be i.i.d. pre and post change but will have some 
dependencies due to ISI (intersymbol interference). However this dependence will be weak because only 
a few symbols at the OFDM symbol boundary will get affected by the symbols of the previous OFDM 
symbol. Thus, one can continue to assume that {V{j,l},j > 1 is an i.i.d. sequence, which is required 
to obtain the simplified algorithm described above. However the i.i.d. may not hold for the CP detector 
because CP resides near the boundary only. Thus, this case will require further consideration. However, 



17 

we will see later, that in the sequential setup, energy detector significantly outperforms the CP detector 
in all possible scenarios we consider. 

The performance of DualCUSUM has been obtained theoretically in [20] and [22]. It is a very 
efficient algorithm because it uses CUSUM at the local cognitive detectors and at the fusion node. 
Also, the local nodes transmit to the fusion node only if they are convinced that there is a change. This 
minimizes cognitive transmissions to the fusion node resulting in low transmit power consumption from 
cognitive nodes and low interference to the primary. Physical layer fusion (see (l22l) ) at fusion node (i.e., 
simultaneous transmissions from all cognitive users) further reduces this interference and also reduces the 
Expected Detection Delay (EDD). Its comparison with several other existing spectrum sensing algorithms 
is available in [14]. 

B. DualCUSUM with Timing Offset 

With an unknown timing offset, the decision statistic used at each node for the CP detector is as follows. 
First, the timing offset estimator of ^ is not preferred here, as under low SNR conditions, to minimize 
the estimation error, we need a large number M of OFDM symbols [18]. This will mean that the amount 
of memory required will be large. Thus, we propose the following. Each node runs Ld CUSUMs for 
each possible timing offset of the primary. In slot j, each node / computes for m G {0, 1, 2, ..., Ld — 1}, 

Rrij, l,m) = Rel^Yl X{{j - 1)L, + i + m, l)X*{{j - 1)L, + i + m + Ld,l)\ , 
{2alRr{j,l,m)-al,) 

l;j,l,m — n 4: /T ' ^j,l,m — {yyj-l,l,m + <ij,l,m) j K^ ' ) 

Wj,i = , max Wj^i,m, yj,i = 61(H/^,>^). (28) 

|mGU,l,..-Ld — 1} 

This algorithm can be intuitively understood as follows. Before change, all the CUSUMs will typically 
be zero as E[Rr{j, 1,771)] = before change. Once the primary arrives, the CUSUM corresponding to 
the correct timing offset m = 9, will start increasing the fastest as it will capture the correct window of 
length Lc- This is similar to the Generalized Likelihood Ratio algorithm discussed in Section V.C for the 
unknown timing offset 9, where 9 ^ {0,1, ... ,Ld — 1}. More comments will follow in Section V.C. 

None of the impairments at the secondary nodes studied above has any effect at the statistics of 
observations at the fusion node. We assume that the cognitive network knows its channel gains and has 
a better control over its system (this is a commonly made assumption in CR). Thus, the DualCUSUM 
at the fusion node remains unchanged. Furthermore, in our implementation, in slot 1 , each node initially 
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captures Lg + L^ samples. From then onwards, each node captures only Ls samples and uses the last 
Ld samples from slot j — 1 to calculate Rr{j, l,m). It can be shown that Rr{j — l,l,m) and Rr{j, I, m) 
remain uncorrelated. This is because a sample in a set of consecutive L^ samples will be correlated with 
some sample in slot j — 1 or slot j, but not both. The performance of this algorithm is provided in Table 

For the energy detector, since timing offset does not affect the decision statistics as discussed in Section 
IV.A, the algorithm remains the same as in Section V.A. Its performance is illustrated in Table IVIIII A 
minor degradation in performance is observed. This is because in a change-detection setup, the presence 
of a timing offset implies that in the slot the primary comes on, the mean energy is less than a"^ + a"^. 

C. GLR-CUSUM with Timing Offset, Frequency Offset and Primary Power Unknown 

Now we assume that a'j ^ is unknown. Additionally, timing and frequency offset could also be present. 
Thus, for the CP detector, we cannot use Rr and need to use R as the decision statistic instead (recall that 
Rr = Real{i?} and Ri = Imag{i?}). It is easy to see that when frequency offset is present, post change, 
Rr ~ A^((T^^ cos(27r(/)),(T^/2Lc) and Ri ~ N{a'^ ism{2iT(/)),a^/2Lc). Also, as we have no knowledge 
of primary signal power o"^; , we now have a composite post change hypothesis, hence we use the 
GeneraUzed Likelihood Ratio (GLR)-CUSUM algorithm [11]. 

The GLR algorithm is briefly described as follows. Let /o be the density of the decision statistic 
Xj^i before change and let fg be the density post change. Here is a parameter that characterizes the 
post-change distribution. In the case of CUSUM algorithm, the parameter 9 is known and the CUSUM 
algorithm in slot j can be described as 

Ty,,= maxfx:iogf^^)V (29) 

Equation (|29l ) can be shown equivalent to (|2TI ). In the case of GLR algorithm, 6 is unknown, but 9 G 
C K, where 3? denotes the real line. Thus ( [29l ) is changed to 

W,,i = max (sup ^ log (^§4) ] • (30) 

In Section V.B, for the unknown timing offset scenario, the algorithm implemented can be described 



as 



T^^l = mf{k : max max ^^ log ( / '''\ ) > t}- (31) 
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It should be noted that here sup is replaced with max as the set is finite and the max over the unknown 
parameter 9 is moved outside. This is done because in the unknown timing offset scenario, keeping the 
max operation inside complicates the computations and requires much larger window sizes for the CP 
detector. This interchange possibly compromises the performance. However, from our simulations we 
will see that the degradation is negligible. 

Now, returning to the current impairments in OFDM, namely unknown frequency offset and primary 
signal power, the supremum is explicit. This is obtained by differentiating the likelihood ratio with respect 
to the unknown (y1i,ip and equating it to zero, and finally substituting the cr'^i,4' which maximizes the 
likelihood ratio. Thus, the GLR test in combination with an unknown timing offset is given by 

Wj^l^rn = max ^ , ,. A /T ' 

i<*<i {j-t + l)ai/Lc 

where Rj. and Ri are the real and imaginary parts of 

1 ^' 
R{j, l,m) = —Y^ X{jLs +i + m, l)X*{jLs +i + m + Ld,l), m £ {0, 1, 2, ..., La - 1}. (32) 

^^ i=i 
The above equation can be intuitively understood as follows. Before the change mean of both Rr and Ri 

are zero, and thus Wj^i^rn will be close to zero. After the change, for the m = 6 (i.e., for the CUSUM 

corresponding to the correct timing offset) since the mean is nonzero, Wj^i^m will keep increasing with 

j, thus eventually detecting the change. The rest of the steps at each secondary node are the same as in 

Section V.A. At the fusion node, the DualCUSUM operation remains unchanged. The computations in 

the GLR algorithm can be limited to a finite window as suggested in [11]. 

For the energy detector, the frequency offset does not affect the performance but due to lack of 

knowledge in primary power we need to use the GLR algorithm. The energy in each slot is 

y^ f^mu^:2)h±M}l^A, (33) 

i=l 

(Here subtraction by o"^ is performed for convenience, for making mean zero before change). At slot, j, 
the GLR algorithm is as follows: 

Wj I = max Aij I, where 

j j 

p=i p=i 1 / CT^ \ /'-! A\ 
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And 6i{i,j,l) is obtained by solving the quadratic equation for 9i 

U-i + l)0l + 01 (2(j - i + l)al + ML,{j - i + l)al + Si,,y) 

- {MLsSQ,^,^i + ML,alS,^j^i - {j - i + l)at) = 0. (35) 

3 3 

where S'j ,, / = N^ V{j^ I) and SQij^i = 2^ ^(j) 0^- Ii^ the above equation 9i{i,j, I) denotes an estimate 

p=i p=i 

of (T^ I (assuming primary has come ON in slot i) and is chosen from 6i G [0, oo). The quadratic equation 
for 9i was obtained by simply differentiating the likelihood ratio w.r.t 9i and setting equal to zero. The 
rest of the steps at a secondary node are same as in DualCUSUM, and fusion node continues to use the 
CUSUM algorithm. The performance of this algorithm is illustrated in Table IVIIII 

D. Algorithms for all Impairments 

We assume that all the above mentioned impairments (including IQ imbalance) could be present and 
(7^ and a1 are unknown to the secondary nodes. For the CP detector, while we can extend the GLR test 
to cover this scenario as well, we have found via simulations, that it is better to first compensate for the 
IQ-imbalance in each slot using (IT2l ) and (IT3]) . Then we estimate the noise power as 

3 L, 

5;5;|X((p-l)L, + fc,0l^ 

(^.n,j,l = -^^ • (36) 

This approximation is valid under low SNR assumptions, as we assume the same value for the variance 
under either hypothesis. Now, since the IQ imbalance can be assumed to have been corrected and we 
have an estimate of noise power ct^ ,■ /, we can use the setup of Section V.C for the other impairments 
(timing offset, frequency offset and received primary power). Thus, for CP detector each node, does the 
same as in Section V.C using the estimated noise power a^ ^ in slot j. 

For the energy detector, the uncertainty in noise power requires the modified GLR (MGLR) algorithm 
[11]. In a CP based detector this was not required as it performed detection of change in the mean of a 
Gaussian signal, and before the change, the mean was known to be zero. Thus, the unknowns are post- 
change mean, and the variances before and after change. Since, under low SNR, variance approximately 
remains same before and after the change, GLR can be used as discussed in the previous paragraph. But 
in case of the energy detector, while the unknown variance, is approximately the same (under low SNR) 
before and after change, the mean both before and after change is also unknown and thus we need to use 
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a modified version of GLR (MGLR) algorithm. To clarify a bit more, in comparison to (BOt the MGLR 
equation will look as 



k 



T-y,i =inf{k : max sup Vlog (/^'(Xi^i)) + sup V log (/e"(^i,0) 

i<s<M>A/- Ve'get^ J \e"eei^^ J 

-(supVlog(/,(Xu))) >7} (37) 

where 6' and 6" are possible parameters before and after change, and 9 evaluates the possibility that there 
is no change. The MGLR approach was first outlined in [11]. The method relies the presence of M* 
samples pre-change. Loosely speaking, the initial set of samples where Hq is the true hypothesis, helps 
in estimating the unknown parameters for subsequent sequential detection of change in the presence of 
impairments. The value of M* depends upon the minimum SNR at which we need to detect reliably and 
the PpA desired. The MGLR algorithm for the energy detector becomes 

^i.' = ,P^^ .^i,jAe^{i+ld.l)>e^il,^,l), whcrc, 



MLs 



Y,iVip,l)-ei{a,b,l)y 



A..i = B\(l) + BU,il)-Bm, i^^(0 = -^%,^(^^,^,).^^^^ • (38) 

And 6i{a, b, I) is obtained by solving the quadratic equation for 6i 

{b-a + l)9l + 9iML,Sa,b,i - MLsSQa,b,i = 0. (39) 

In the equation ( [39l ) 6i{l,i,l) is an estimate of a'^ ^ and 6i{i + l,j,l) is an estimate of a'^^ + "^^ ;> 
assuming primary has come on at slot i + 1. 9i{l,j,l) is an estimate of a'^ assuming primary has not 
come on. The rest of the steps at a secondary node are the same as in DualCUSUM and fusion node 
continues to use CUSUM. The performance of this algorithm is illustrated in Table IVIIII 

The condition in ( [38l ) is for detecting OFF— )-0N, i.e., we are detecting an increase in signal power. 
The condition needs to be reversed for detecting ON— >OFF [11]. 

E. Performance Comparison 

In this subsection, we compare the performance of the above algorithms. There are 5 nodes. The SNR 
at each node is — lOdB. Wherever applicable, (j) = 0.1, Ac/) = 10°, e = 0.2,6 = 10. The change time T 



22 



Pfa 


(IV.A) 


(IV.B) 


(IV.C) 


(IV.D) 
(all impairments) 


Snapshot 
(all impairments) 


0.1 


10.15 


18.27 


24.71 


28.15 


64.16 


0.075 


11.43 


19.82 


28.07 


31.01 


67.46 


0.05 


12.6 


22.09 


31.42 


34.95 


72.35 



TABLE VII 

CP BASED CO-OPERATIVE SPECTRUM SENSING ALGORITHMS. 



Pfa 


(IV.A) 


(IV.B) 


(IV.C) 


(IV.D) 
(all impairments) 


Snapshot 
(all impairments) 


0.1 


5.22 


5.43 


7.73 


10.15 


349.13 


0.075 


5.61 


5.91 


8.52 


11.43 


438.03 


0.05 


6.41 


6.46 


9.19 


12.6 


623.58 



TABLE VIII 

Energy detector based co-operative spectrum sensing algorithms. 



(in units of OFDM symbols) is assumed to have a geometric distribution with parameter p = 0.004. For 
different values of Pfa, EDD in units of OFDM symbols is shown in Table IVlIl for CP-based detectors 
and in Table I VIII I for energy-detector based algorithms. 

For comparison, we have also simulated a cooperative snapshot detector for both CP and energy 
detectors. CP detector captures M = 50 OFDM symbols of data and detects the signal in the presence 
of all impairments and compensating for the same using the steps in Section III.E. The energy detector 
captures M = 5 OFDM symbols and detects the signals in presence of all impairments. Compensation 
is done for IQ imbalance. The values of M chosen for the two snapshot detectors are chose to minimize 
EDD in each case for a given Pfa- Each node sends a 1 or according to whether Hi or Hq is chosen. 
The fusion center uses the AND rule to decide between Hq or Hi as the AND rule works the best in 
the present setup. For the snapshot detector, we assume that the fusion node has no noise. 

We see that, as the amount of uncertainty increases, the performance degrades for both the detectors. 
Also, from last two columns, we see that the sequential setup provides significant performance gains over 
the snapshot detector (even though for the snapshot detector we have assumed no noise at the fusion 
node) for both the CP-detector and the Energy-detector. Also, comparing the energy detector and the CP- 
detector we can clearly see that in a change-detection setup, energy detector significantly outperforms 
the CP detector (by comparing the columns labeled IV.D in both the tables) under all scenarios. (In this 
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example M* for the MGLR was chosen as 50 OFDM symbols). But the snapshot energy detector shows 
significant degradation under noise uncertainty. 

VI. Conclusions 

We have considered the problem of spectrum sensing of OFDM signals using cyclic prefix based 
and energy based detectors. We have analyzed the effect of some typical impairments like timing and 
frequency offset, IQ-imbalance and transmit/noise power uncertainty and presented techniques to modify 
the detectors to work under these impairments. We have also proposed cooperative sequential change 
detection based algorithms and overcome the effects of these impairments in that setup also. We have 
shown that sequential detection improves the performance significantly as against fixed sample size 
detectors. It is also shown that the sequential energy detector significantly outperforms the CP detector 
under all impairments but the snapshot energy detector performs worse than the CP detector under noise 
power uncertainties. 

Most of these detectors will work under time varying multipath frequency selective fading also. In 
future work we will verify these claims via simulation and also consider frequency selective fading under 
different impairments discussed in this paper. 
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